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Abstract 

We have recently reported the calculation of slope and curvature 
of Isgur-Wise function based on Variationally Improved Perturbation 
Theory (VIPT) in a QCD inspired potential model. In that work, 
Coulombic potential was taken as the parent while the linear one as 
the perturbation. In this work, we choose the linear one as the parent 
with Coulombic one as the perturbation and see the consequences. 
Keywords: VIPT, Isgur-Wise function, charge radii and convexity pa- 
rameter. 

PACS Nos. 12.39.-x ; 12.39.Jh ; 12.39.Pn 



1 Introduction 

The use of VIPT [1] in the calculation of slope and curvature of Isgur-Wise(I- 
W) function is very recent. In a QCD inspired potential model , we have 
two options in using VIPT [21 |3l H] with the linear cum Coulomb poten- 
tial of QCD[18] -(i) Coulombic one as the parent, the linear one as the 
perturbation and (ii) linear one as the parent, Coulombic one as the per- 
turbation. Already, Coulombic part as parent has been anaiyzedpQ succesfully 
taking into account of three terms in the summation of equation(lO) of Ref 
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1 for D, D s andB mesons .However, the results were shown to be improved 
only with more terms in the equation(lO) of Ref l.In this work, we take an 
alternate strategy : we take the linear potential as the parent and see any 
modification over the results of earlier workjT] for the same three number of 
terms of equation(10)in that work. We recall that [3] for the linear potential 
to be dominant we prefer < r > > r ,where < r > is the expectation value 
of the distance r which reasonably gives the size of a state (in this case me- 
son) and ro is a point at which linear cum Coulomb potential becomes zero( 
Fig.l of Aitchison and Dudek ,Ref [3] ).In this work we have also included 
the B s andB c mesons which have greater reduced mass to check the range 
of applicability of the approach. 

The rest of the paper is organized as follows : section 2 contains the 
formalism, section 3 the result and calculation while section 4 includes the 
discussion and conclusion. 



2 Formalism 

2.1 Isgur-Wise function; its slope and curvature 

The Isgur-Wise function is written as [5] : 

= ay) 

= l-p*(y-l) + C(y-l) 2 + ... (1) 

where 

y = v< (2) 

and and v' being the four velocity of the heavy meson before and after 
the decay. The quantity p 2 is the slope of I-W function at y = 1 and known 
as charge radius : 

The second order derivative is the curvature of the I-W function known as 
convexity parameter : 



For the heavy-light flavor mesons the I-W function can also be written as 

[sum : 

£ (y) = / 4irr \ij)(r)\ cos prdr (5) 

J 

where 

p 2 = 2/i(y-l) (6) 

Now the wavefunction ip of the hadronic system is determined by taking the 
linear as parent potential. 



2.2 First order corrected wavefunction and energy in 
VIPT 

The wavefunction corrected upto the first order of j th state is given by (Equa- 
tion 10 of Ref 1) : 



The energy corrected upto first order for the same state is 
E 3 = Uf^H^dv 



Jijf*(H oP , + H' pl )ijfdv (8) 



where ipk, are the wavefunction and energy eigen values of the k states 
which are orthonormal to j th state .The superscript (0) means zeroeth order 
correction of the corresponding quantities .Also , we note that P'is the vari- 
ational parameter and H p> , H' P , defined in equation(9)of Ref 1. 

The summation in equation (7) above [Equation 10 of Ref 1] , can include 
any number of k states .In this work , we consider terms upto three states 
in the summation as was done Ref[l]. 
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2.3 Wavefunctions using VIPT with linear potential 
as the parent 

(i)With one term in the summation 

As explained earlier , we take V as the variational parameter instead of the 
physical parameter b in the parent linear potential to write the Hamiltonian 
as PHI?]: 

H = H + H' 

= h or h c 

2fi 3r 

V 2 , a 

= V or he 

2/i r 

= V or or + or + c 

2/i r 

= H ob , + H' b , (9) 

where a = Now ,H Q y = — y- — b'r is the parent Hamiltonian with the 
new parameter 6' and if^' = - — &V + &r + c is the perturbed Hamiltonian 
with the same variational parameter b' instead of the physical parameter b . 

We consider j th as Is state (n = 1,1 = 0) and in the summation of equa- 
tion^) [ equation(lO) of Ref 1], we consider a single k th st&te which is the 2s 
state (n = 2,1 = 0). 

We note that in the variational method, we are interested only in the V 
dependence of the Hamiltonian, and so ' d in H' b , has no role to play in the 
calculation [9]. 

The unperturbed wavefunctions with linear parent with appropriate bound- 
ary conditions are the Airy functions given by [3]: 

VV) (r) = ^~^Ai f(2/x6')^ + Pan) (10) 



20rr 

where po n s are the zeroes of the Airy function Ai (p on ) = given by [6]: 

[ 37r(4n-l) ^ 

Pon = ~ Z l 11 ) 
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Table 1: A few of the zeroes of Airy function for s states . 



State 


POn 


Is 


-2.3194 


2s 


-4.083 


3s 


-5.5183 


4s 


-6.782 



and N n is the normalization constant. 

As an illustration , we reproduce for s states a few of the zeroes of the 
Airy function in table 1. 

The corresponding energies are given as : 

£„=-(y%„„ (i2) 

Of course n — 1, 2, 3, 4, is the principal quantum number. 

Thus the trial Is state (n — 1,1 — 0) wavefunction is (which is also the 
unperturbed wavefunction) : 



(0) 
10 

r 

Ai( Zl ) (13) 



= — ^Ai((2ub')*r -2.3194") 



2^Txr 



where 



^ = ((2^)^-2.3194) (14) 
and the subscript 10 indicates the quantum number [n, I) of the j th state. 

We note that b' is replaced by b' which is obtained by minimizing Ej 
given by equation(8).It is essential since in VIPT we have to use the values 
of variational parameter leading to minimum energy (for example in Ref 1, a s 
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was replaced by a' 10 ). The values of b for different mesons are listed in table 2. 

Now we consider the single k th state in the summation of equation (7)which 
is the 2s state given by : 

4? = ^M^ir- 4.083) 
No 

' At(z 2 ) (15) 



2^/ir 
where 



Z2 



(2//6')3r- 4.083) (16) 



The wavefunction corrected upto first order is : 



ij s = N 



^ + ~r ((& - 1) < r > 2>1 -a < i > 2>1 ) ^ 20 (r) 

(po2-poi)b 

(17) 

where 

< r > 2 ,i= A^iiV 2 ^ + °° ((2/jb')^r - 2.3194) ((2^)3r - 4.083) (18) 
and iV is the normalization constant. 



(ii)With two terms in the summation 

We next consider the 3s state (n = 3,1 = 0) in addition to 2s state (as done 
in the single term case) given by : 

4? = ^^((2^)^-5.5153) 



2, 

N 3 

where 



Ai(z 3 ) (19) 



z 3 = ((2/x6')3r- 5.5153) (20) 
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With the inclusion of this state , the wavefunction corrected upto the first 
order is : 

ipD = N'[^ + (2/i)3 ((& -b')<r > 2)1 -a < l - > 2)1 ) ^ 2 o (r) + 

(P02 - Poi) b 

1 

(2 ^ )3 _ /§ ((6 - < r > 3>1 -a < ± > 3;1 ) ^30 (r)](21) 

(P03 - P0l) & 

where 

< r > 3i i= A^iiV 3 ^ + °° Ai ((2/i^)3r - 2.3194) ((2^)^ - 5.5153) (22) 
and A/"' is the normalization constant. 



(iii)With three terms in the summation 

In addition to the 2s and 3s states we now add the 4s state : 

4? = -^^((2^)^- 6.782) 

= < 23 » 

where 

z 4 = ((2fiV)sr- 6.782) (24) 
With the inclusion of this state , the first order wavefunction now becomes : 

i 

ipT = N"[^ + (2/i)3 _ A ((& -b')<r > 2>1 -a < I > 2>1 ) ^20 (r) + 

(P02 - Poi) b 



(2p) 



(P03 - P0l) 

(2/i) 



-|- ((& - b') < r > 3 ,i -a < J > 3 ,i) ^30 (r) 



(Po4 - An) & 



-j (jb - 6') < r > 4 ,i -a < J > 4 ,i) ^40 (r)](25) 



where 

< r > 4 ,i= N^J^ Ai([2^)^r- 2.3194) Ai{{2^)h- 6.782) (26) 

and N" is the normalization constant. 

The relativistic version of these wavefunctions are obtained by multiplying 
above expression by (rfia)~~ e [HI [12]. Thus, relativistic version of all these 
wavefunctions is: 

A,rd = ipi {rna)~ e (27) 

where i — S,D,T 
and 

Putting all these wavefunctions i.e. equations (17), (21), (25) and (27) in 
(5) we can calculate the Isgur-Wise function for the different cases. 



3 Calculation and Results 

We have listed the values of charge radius and convexity parameter of the 
calculated I-W function for various heavy-light flavor mesons in the present 
method considering single state , two states, and three states in the summa- 
tion occurred in VIPT with perturbative Coulombic and relativistic effect. 

In tables(3-5), we record the predictions of the present work for single 
term,two terms and three terms respectively .In table 6, we give comparision 
of this work to that of Ref[l]; while table 7 gives a summary of the predic- 
tions of slope and curvature of IW function in other models. 

The a s values are taken from the l^-scheme [331 El HHJ EE] and the inte- 
grations are done numerically for all these calculations. 
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Table 2: Values of b' with b = 0.183GeV 2 . 



Mesons 


Reduced mass /i 




b' without relativistic effect 


fr'with relativistic effect 


D 


0.2761 


0.924 


5.406 


16.24 


D s 


0.368248 


0.924 


5.876 


19.8 


B 


0.31464 


0.348 


4.33 


5.587 


B s 


0.4401 


0.348 


4.497 


5.954 


B c 


1.1803 


0.348 


5.39 


8.103 



Table 3: Values of slope p 2 and curvature C with single term in equation(7) 



Meson 


P 2 s 


C s 


Ps,rel 


Cs,rel 


D 


1.36 


0.01 


0.53 


0.0022 


D s 


1.867 


0.03 


0.702 


0.0036 


B 


1.93 


0.02 


1.41 


0.013 


B s 


2.923 


0.046 


2.113 


0.0283 


B c 


9.442 


0.484 


6.274 


0.2522 



Table 4: Values of slope p 2 and curvature C with two terms in equation(7) 



Meson 


Pi 


c D 


p 2 D,rel 


Co,rel 


D 


1.201 


0.013 


0.57 


0.0026 


D s 


2.001 


0.0242 


0.74 


0.0041 


B 


2.004 


0.0244 


1.44 


0.0133 


B s 


3.031 


0.0565 


2.16 


0.0297 


B c 


10.2 


0.61 


6.51 


0.275 
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Table 5: Values of slope p 2 and curvature C with three terms in equation(7) 



Meson 


Pt 


Cj 1 


PT,rel 


Cr,rel 


D 


1.33 


0.016 


0.604 


0.00326 


D s 


2.023 


0.0305 


0.78 


0.0054 


B 


2.027 


0.031 


1.54 


0.0217 


B s 


3.087 


0.071 


2.29 


0.047 


B c 


10.25 


0.767 


6.99 


0.441 



Table 6: Comparision of this work to that of Ref[l] with relativistic effect 
being included. 



Me- 
son 


Ref[l] 


This work 


Ps 


C s 


Pi 




Pt 


Ct 


Ps 


C s 


Pi 


c D 


Pt 




D 


.433 


.525 


.432 


.524 


.43 


.516 


0.53 


0.0022 


0.57 


0.0026 


0.604 


0.0033 


B 


.56 


.85 


.55 


.84 


.545 


.815 


0.702 


0.0036 


0.74 


0.0041 


0.78 


0.0054 


D s 


3.6 


15.3 


3.16 


12.32 


3.12 


11.8 


1.41 


0.0126 


1.44 


0.0133 


1.54 


0.0213 


B s 














2.11 


0.0283 


2.16 


0.0297 


2.29 


0.0471 


B c 














6.27 


0.252 


6.51 


0.275 


6.99 


0.441 
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Table 7: Predictions of the slope and curvature of the IW function in various 
models. 



Model 


Value of p 2 


Value of curvature C 


Le Youanc et al [19] 


> 0.75 




Le Youanc et al [20] 


> 0.75 


> 0.47 


Rosner [28] 


1.66 


2.76 


Mannel [291 ED] 


0.98 


0.98 


Pole Ansatz [31] 


1.42 


2.71 


MIT Bag Model [27] 


2.35 


3.95 


Simple Quark Model [26J 


1 


1.11 


Skryme Model [24] 


1.3 


0.85 


QCD Sum Rule [25] 


0.65 


0.47 


Relativistic Three Quark Model [23J 


1.35 


1.75 


Infinite Momentum Frame Quark Model [22J 


3.04 


6.81 


Neubert [32] 


0.82±0.09 





4 Discussion and Conclusion 

We note that in the present analysis , the slope and curvature increase with 
the inclusion of more terms in the summation of equation (7) in contrary to 
the case of Coulombic parent [1], where with more terms the slope and cur- 
vature decreased. An analysis of table 6 also indicate that for a definite term 
, the slope has increased in the present work compared to Ref[l], while for 
the curvature ,the pattern is reversed i.e. it has decreased in the present 
case compared to Ref[l].So,with only one term for the linear parent in equa- 
tion^), we can have comparable values for slope and curvature with other 
models and data [table- 7]. This is undoubtably a great phenomenological ad- 
vantage as involvement of more terms in equation(7) makes the calculation 
quite cumbersome which happened in Ref[l]. However, relativistic correction 
in this case also decrease the slope and curvature of Isgur-Wise function as 
observed earlier [TJ. 

To conclude , the present approach with the linear potential as the par- 
ent has provided much more improved results for the slope and specially 
curvature of the IW function within the QCD inspired potential model mi- 
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der study and appeared to be preferable over the one of Ref[l] where the 
linear potential was considered as perturbation. 
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